HD-TVP-94-16 



How sharp is the chiral crossover phenomenon for 

realistic meson masses? 

Hildegard Meyer-Ortmanns* 
^ ! and 
C\ ' Bernd-Jochen Schaefer^ 

O-i- 

<D ; Institut fiir Theoret. Physik, 

Universitat Heidelberg 
(N Philosophenweg 19 

D-69120 Heidelberg, F.R.G. 

O 

m 

o\ 
o 

' The mass dependence of the chiral phase transition is studied in the linear 

^ \ SU(3) x SU(3) sigma-model to leading order in a 1/iVj-expansion, Nj de- 

noting the number of flavours. For realistic meson masses we find a smooth 

Ph. 

crossover between T ~ 181.5 to 192.6 [MeV]. The crossover looks more rapid 



X 



Abstract 



in the light quark condensate than in thermodynamic quantities like the en- 
ergy and entropy densities. The change in the light quark condensate in this 
temperature interval is ~ 50% of the zero-temperature condensate value, while 
the entropy density increases by (5.5 ±0.8) • 10 -3 [GeV 3 ]. Since the numerical 
error is particularly large in this region, we cannot rule out a finite latent heat 
smaller than 0.2 [GeV/fm 3 ]. The chiral transition is washed out for an average 
pseudoscalar meson octet mass of 203 [MeV]. This gives an upper bound on 
the first-order transition region in the meson mass parameter space. The cor- 
responding ratio of critical to realistic light current quark masses mffj j m Ut d 
is estimated as 0.26 ± 0.08. This result is by an order of magnitude larger 
than the corresponding mean-field value. Therefore the realistic quark/meson 
masses seem to lie less deeply in the crossover region than it is suggested by 
a mean-field calculation. 
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1 Introduction 



Spontaneous symmetry breaking is frequently used as an ansatz to explain the driv- 
ing force of finite temperature phase transitions in QCD. The symmetries refer to 
certain limiting cases of the QCD-Lagrangian. In the limit of infinite quark masses 
and in the case of three colours finite temperature QCD is invariant under Z(3)- 
transformations (Z(3) is the center of SU(3)). The spontaneous breaking of Z(3) is 
made responsible for the phase transition from the confinement phase at low temper- 
atures to the deconfinement phase at high temperatures. In the other extreme case 
of Nf vanishing current quark masses QCD is invariant under SU(Nf) x SU(Nf) 
chiral transformations. The restoration of the spontaneously broken chiral symme- 
try at high temperatures is said to drive the chiral phase transition of QCD. For the 
chiral limit the renormalization group analysis of Pisarski and Wilczek || serves as 
a guideline for conjectures about the order of the chiral phase transition. Similar 
studies have been performed by Svetitsky and Yaffe || in the other extreme case of 
pure gauge theory or infinite quark masses. 

In reality the current quark masses are neither infinite nor zero. The masses of 
the charm, bottom and top quarks are large compared to the scale of the critical 
temperature T c of QCD, which lies between 150 and 250 [MeV]. Thus it seems to 
be justified to treat them as infinite in thermodynamic investigations. The (renor- 
malization group invariant current) masses of the up and down quarks are small 
compared to T c , m u = 7.6 ± 2.2 [MeV], m d = 13.3 ± 3.9 [MeV] for A = 100 [MeV], 
where A refers to the MS scheme with 3 flavours ||. At a first glance it seems to be 
well justified to set them simply to zero, although more careful investigations in the 
framework of chiral perturbation theory give us a warning to neglect their influence 
on the chiral transition [|J. A particular role is played by the strange quark. Its 
mass (m s ~ 205 ± 50 [MeV]) is neither small nor large, but just of the order of the 
scale which is set by T c . To get an idea about the influence of finite quark masses 
on the phase structure of QCD it is instructive to consult statistical physics. The 
phase transition from a liquid to a gas is of first order below some critical value of 
the pressure. For a critical value of the pressure it becomes of second order. Above 
the critical strength it turns into a smooth crossover between the liquid and the 
gas phase. Similar effects are known from ferromagnets under the influence of an 
external magnetic field. In an 0(A r )-ferromagnet an arbitrarily small magnetic field 
is sufficient to turn the second order transition with an infinite correlation length 
into a crossover without a diverging correlation length. In close analogy to our 
subsequent considerations we should also mention the three-dimensional Z(3)-Potts 
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model. The spin variables of this model can take three values at each lattice site. 
For a vanishing external magnetic field the Potts model is known to have a first 
order phase transition at finite temperature. For a critical field strength it becomes 
of second order and disappears for even stronger external fields. 

The conjugate variables are the pressure and the specific volume in the liquid/gas 
system, the external magnetic field and the magnetization for the ferro magnet. The 
analogous pair in QCD are the current quark masses and the quark condensates as 
the associated order parameters for the chiral transition (on the quark level). On 
the mesonic level we have external fields, which can be related to the quark masses, 
and mesonic condensates. The comparison between the statistical systems and QCD 
goes beyond a formal analogy. In an SU (3)-lattice gauge theory it has been shown 
by DeGrand and DeTar || and Banks and Ukawa || that dynamical fermions on 
the QCD level induce an external field coupled linearly to a spin field in an effective 
Z(3)-spin model. The spin model has been derived from QCD in a strong coupling 
expansion at high temperatures. The external field strength was expressed in terms 
of the hopping parameter of the lattice-QCD Lagrangian. 

In this paper, we consider the SU (3) x SU (3) linear sigma-model as an effective 
model for chiral symmetry restoration. We include two external fields to account 
for the finite quark masses on the mesonic level. The action is constructed in terms 
of the chiral order parameter field, which is a matrix of mesonic condensates. In 
the special case of one external field and one mesonic condensate, the action of 
the sigma-model takes the form of Landau's free energy functional for a liquid/gas 
system (cf. section 5). Thus we expect the same qualitative features as we know 
from the liquid/gas system. Beyond certain critical values of the meson masses the 
chiral transition should be washed out and turn into a crossover phenomenon with 
smooth changes in the condensates, energy density and entropy density. Therefore 
the main question is, as to whether the realistic quark masses are too large for the 
chiral transition and too small for the deconfinement transition to persist. 

Numerous numerical simulations in lattice QCD have investigated the phase 
structure of QCD. Results on the mass dependence of the finite temperature transi- 
tions were summarized in the famous Columbia-plot J7j, cf. Fig. 1. Indicated are the 
presence or absence of QCD transitions as a function of the bare quark masses m s -a 
and m Uj( i ■ a in lattice units, a denotes the lattice constant. The solid line stands for 
second order transitions, the shaded areas enclose mass values leading to first order 
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Figure 1: Columbia plot [7]. Generic phase diagram, partly conjectural, for 2 and 3 
flavours. Solid lines indicate supposed 2nd order transitions, shaded areas 1st order 
transitions. Solid circles correspond to mass parameters, which lead to 1st order 
transitions, triangles to mass parameters, which lead to crossover phenomena. The 
open circle locates the suggested physical mass point. 

transitions. Solid circles correspond to mass values leading to 1st order transitions, 
e.g. the point with 3 degenerate bare quark masses m u a = m^a = m s a = 0.025. 
Triangles indicate mass combinations which lead to a crossover phenomenon. In 
particular, the mass point is included with two light (m u a = m^a = 0.025) and one 
heavier flavour (m s a = 0.1). This mass point comes closest to the realistic quark 
masses with a ratio of m s /m u ^ ~ 20. The ratio between the critical quark mass 
m crit which ii es on the first order transition boundary to the realistic quark mass 
m u is about 0.5 if the ratio m^a/mf^a = m u a/m s a is kept fixed. Fig. 1 suggests 
that the masspoint with realistic quark masses is not too far from the second or- 
der transition line, although the precise values of the critical strange quark masses, 
where the first order transitions turn into second order transitions, are still an open 
question. 

Thus the conclusion from the lattice is that there is no true chiral phase transition 
for experimental quark masses. Depending on the degree of smoothness of the chiral 
transition this conclusion can have far reaching consequences in view of measurable 
effects in heavy-ion collisions. Therefore we spend some comments on the reliability 
of the lattice result. 

Monte Carlo results from lattice calculations are in general affected by artifacts 
due to the finite IR-cutoff (i.e. the finite volume), the UV-cutoff (i.e. the finite 
lattice constant) and finite bare quark masses (it is unavoidable to extrapolate to 
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the chiral limit). Applied to the masspoint with two light and one heavier flavour 
{rriu,dCL = 0.025, m s a = 0.1) of the Columbia plot [7j there are some warnings which 
should be kept in mind. The lattice extension in imaginary time direction was 
chosen as N T = 4. Even without dynamical fermions N T should be larger than 10 
to find a critical bare coupling close to the continuum limit. One manifestation of 
the vicinity to the strong coupling regime are the masses of flavour partners which 
should be degenerate in the continuum limit. The ratio of two .fT-meson masses is 
still about 2 |7]]. The effect of heavier flavour partners on the order of the chiral 
transition is difficult to control. 

Another subtle effect of the finite UV-cutoff is a change in the effective symmetry 
as a function of the bare coupling if staggered fermions are used |§. For four 
continuum flavours the chiral transition is of 2nd order in the strong coupling limit, 
but of 1st order close to the continuum region. For two flavours the critical indices 
of the second order transition change between the strong and weak coupling region 
|J. For three flavours one can speculate, whether the crossover itself is a strong 
coupling feature rather than continuum physics. 

A further obstacle in the staggered fermion formulation applies in particular to 
the case of three flavours. The effective fermionic action which projects on three 
flavours must be regarded as a prescription rather than being derived from the 
staggered fermion action J7|. A derivation could only lead to integer multiples of 
four continuum flavours in the effective action. The error due to this prescription 
seems to be difficult to control. 

A bulk "chiral" transition has been identified for 8 flavours |TI| (for bulk transi- 
tions the transition itself is a mere lattice artifact). Precursors of this bulk transition 
for a smaller number of flavours can be superimposed on the finite T-chiral transition 
and influence its strength. 

If the bare coupling lies in the strong coupling regime, the conversion of bare 
masses in lattice units into masses in [MeV] is not unique. Last but not least it is not 
clear whether the ratio of m u ^a/m s a is unaffected by multiplicative renormalizations 
(it may be that / m r s en ^ m u ^a/m s a). Hence the very location of the physical 
mass point in the (m Ut da, m s a) diagram is questionable. 

These warnings should suffice as arguments in favour of alternative approaches which 
include different approximations to study the phase structure of QCD. 



In [11] a systematic study of the mass dependence of the chiral transition has 
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been proposed in the framework of effective models for QCD. In a first paper [|Tl|] a 
bound on the first order transition region was given for the 5'?7(3)-symmetric case. 
The loose bound on the (average) pseudoscalar octet mass was 100 [MeV] for a 
sigma meson mass between 600 and 950 [MeV]. The qualitative conclusion was that 
realistic meson masses lie deeply in the crossover region. The results were obtained 
in a large iV^-expansion under the omission of n ^ 0-Matsubara frequencies. 

In this paper we will give a more stringent bound on the maximal pseudoscalar 
octet mass, for which the transition is still of first order. The large iVy-expansion 
seems to be a good starting point for the SU (3) x SU(3) linear sigma-model, which 
reduces to an 0(18)-model in a special case of two vanishing couplings, see section 
3 below. Therefore we still use the large iV/-expansion, but include also the n^O- 
Matsubara frequencies. The omission of the n ^ 0-modes seems to be justified only 
for high temperatures or in case of a second order phase transition. The distance 
(in mass parameter space) between experimental meson masses and critical meson 
masses, for which the first order chiral transition turns into a crossover phenomenon, 
is more than a quantity of academic interest. Its physical relevance can be seen as 
follows. Assume that the physical masses are not identical but very close to the 
critical mass values where the chiral transition is of second order. In such a case the 
chances are good to see some remnant in heavy-ion collisions of a diverging correla- 
tion length for critical mass values, i.e. a rather large correlation length for realistic 
masses. 

One manifestation of a large correlation length has been supposed to be large clusters 
of charged or neutral pions which are aligned in isospin space [TT2]. The correlation 



volume of a cluster should be large enough that the number of emitted pions is 
sufficient for the detector to resolve the cluster structure. Intermittent behaviour 



has been proposed as another manifestation of a second order transition Hl3] . In 
this case one would observe some kind of self-similarity phenomenon if the bin size 
of rapidity intervals is made finer and finer. 

Also the vicinity to the first order transition region could leave observable effects if 
the crossover for realistic masses is still sharp enough. A sharp crossover is associated 
with a rapid change in the condensates and their induced masses, and/or a rapid 
change of entropy in a small temperature interval. Average transverse momentum 
distributions of charged particles could be flattened as a function of the multiplicity 



dN/dy of final state particles in a given rapidity interval [14|. Pronounced fluctua- 
tions in the particle multiplicities would show up if the crossover is strong enough 
to induce deflagration processes during the phase conversion [[HJ . True singularities 
of a first or second order transition will be anyway rounded in real experiments due 
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to the finite volume. 

It may be academic to ask as to whether chiral symmetry restoration in the infi- 
nite volume limit proceeds via a weakly first order transition or a smooth crossover 
phenomenon. It is certainly more sensible to pose the question in the following way: 
Is the gap in entropy densities in the transition/crossover region sufficient to induce 
multiplicity fluctuations in the observed pion yield, lying clearly above the statistical 
noise? In this paper, we therefore try to find a quantitative measure for how sharp 
the chiral crossover phenomenon is for realistic meson masses. 

The paper is organized as follows. In section 2 we present the tree level parametriza- 
tion of the SU(3) x SU (3)-linear sigma-model at zero temperature. We give a pre- 
scription to translate the meson condensates and meson masses to quark condensates 
and current quark masses. In section 3 we summarize the essence of our approach. 
The mesonic self-interaction is treated to leading order in an expansion in the num- 
ber of quark flavours Nf. On the mesonic level we have 9 scalar plus 9 pseudoscalar 
mesons. The 577(3) x S77(3) sigma-model reduces to an (9(18)-model in a certain 
limit. Thus the leading order of a 1/iV-expansion should be a good starting point. 
The thermodynamic effective potential is evaluated in a high-temperature expan- 
sion and in a fully numerical approach. The numerical approach is certainly more 
appropriate to the phase transition region. To discuss the mass dependence of the 
order of the chiral transition we distinguish three regions in mass parameter space 
(m n , m,K, ■ ■ •) and (m Ut d,m s ): the chiral limit (section 4), several mass points on 
the first order transition boundary, so-called critical mass points (section 5) and 
realistic meson masses (section 6), which come close to the experimental masses. In 
the chiral limit we fix the couplings of the tree level parametrization. The quanti- 
tative discrepancy between the high-temperature expansion and the fully numerical 
evaluation will be revealed. We calculate the barrier height between the coexisting 
minima of the effective potential as one measure for the strength of the first order 
transition. 

Values of critical meson masses are determined in section 5 for three special cases: 
the S77(3)-symmetric case with finite, but degenerate pseudoscalar octet masses, 
the case with a realistic mass splitting between the masses of the pseudoscalar octet 
(the ratios of meson masses with and without strangeness are kept fixed to their 
experimental values), and the case of vanishing strange quark mass m s . The critical 
meson masses are calculated both in a mean-field approximation and in the large 
^-expansion. 

Section 6 deals with meson masses which are rather close to the experimental values. 
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We describe the crossover phenomenon in the meson and quark condensates as a 
function of temperature. The crossover is also manifest in thermodynamic quantities 
like the energy density e, the entropy density s, and the pressure p. We derive e, s 
and p from the partition function of the sigma-model in a saddle-point approxima- 
tion. Upper bounds on a finite latent heat during the chiral transition are predicted. 
The temperature dependence of effective masses which enter the calculation of the 
effective potential will be displayed. In section 7 we summarize our results and draw 
some conclusions in view of heavy-ion experiments. 



2 Tree level parametrization of the SU(3) x SU(3) 
linear sigma-model 

The Euclidean Lagrangian density of the SU (3) x SU (3)-linear sigma-model is given 



as 

Ceuci = ^Tr(d li Md fM M + )-^i4TrMM + + g(detM + detM+) + 

+ A(TrMM + ) 2 + / 2 Tr(MM + ) 2 -5o^o-^8, (1) 
where the (3 x 3)-matrix field M(x) is written as 

1 8 

M=-=^2(<T t +iir t )\t. (2) 

Here ai and lie denote the nonets of scalar and pseudoscalar mesons, respectively, 
\e (£ = 1,...,8) are the Gell-Mann matrices, Ao = J\ • diag(l, 1, 1). The chiral 
SU(3) x S'?7(3)-symmetry is explicitly broken by the term (— eocr — Es&s), which 
is linear in the external fields £o,£s- A non- vanishing value of 6q gives a common 
mass value to the octet of pseudoscalar Goldstone bosons rn n ,mk,rn v . When also 
e 8 7^ 0, it can be adjusted such that it leads to a realistic mass splitting inside the 
(pseudo) scalar meson octets. 

The Lagrangian (|IJ) appears as a natural candidate for an effective model, which 
is designed to describe the phenomenon of chiral symmetry restoration. The action 
S = J d 3 xdrC(x) with C of Eq. fll]) may be regarded as an effective action for QCD, 
constructed in terms of a chiral order parameter field M for the chiral transition. 
As order parameters for the chiral transition we choose the meson condensates (o"o) 
and (as)- The expectation value of M is then parametrized in terms of (ao), (as) 
according to 



(M) = diagi= 



(^o) + -j= (<r 8 ), (a ) + -^(o-s), (ffo) - V2(a 8 ) 



(3) 
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The construction of an action 5 in terms of an order parameter field is a concept in 
close analogy to Landau's free energy functional T in terms of an order parameter 
field, T coincides with 5 in the mean-field approximation. Quartic terms in M have 
to be introduced in Eq. (P to allow for the possibility of spontaneous symmetry 
breaking. For g = = Eq = the Lagrangian is still invariant under U (3) x U(3) 
transformations. For Nf > 3 there are two independent quartic terms, parametrized 
with coefficients fi and f 2 - To account for a realistic (rj, ?/)-mass splitting, a det- 
term with an "instanton" -coupling g has to be included as well. It reduces the 
symmetry of £ to SU(3) x SU(3) if e$ = = e$- Finally, the external field terms 
which are linear in M are the most simple choice for an explicit symmetry breaking 
accounting for the small, but finite masses of the Goldstone octet (m n , tuk , Tn^) . 
Thus one arrives at the 577(3) x SU (3) linear sigma-model in a natural way if one 
is interested in the limited aspect of chiral symmetry restoration. 

It remains to fix the Lagrangian parameters /xjj, fx, f'2, g, £0, £s from an experi- 
mental input. The choice of the input and the way of parametrization are in no 



way unique, cf. fll£], [L7|, [ljj. Since the pseudoscalar meson masses are experimen- 
tally well known, m^, rrik, m v , rn v i and f n have been used in Jll| to fix /xjj, fi, f 2 , g, e 



and e$. In addition, the mass of the cr'-meson has been treated as input parameter 
and varied between 600 and 950 [MeV] to solve for and fi, which occur in the 
(pseudo) scalar meson masses in the combination (— + 4/i(<Tq + erf)). It is worth 
noting that the order parameters (cro), (erg) can be determined without knowing /!q 
and fi separately. The equations for /Xq, fi, f 2 , g, e , e 8 do not admit solutions for 
an arbitrary choice of ■m 7T ,mk,m v ,m rt i , f n and m a . This was the reason why the 
input masses which have been actually used as input in |TT[ were slightly deviating 



from the experimental values if m a , was chosen as 600 or 950 [MeV]. In [T3|] the 
experimental values could be used for m^, m k , m v , m v ', f n on the price that m a , was 
used as input with 1400 [MeV] . 

In this paper, our interest goes beyond the point with (almost) experimen- 
tal pseudoscalar meson masses. Naturally, the question arises why we are not 
satisfied with a parametrization that reproduces the experimental values for the 
(pseudo) scalar meson masses, but want to tune some parameters to unphysical, un- 
realistic mass values as well. There are two reasons for that. The first one is to 
check how stable our results about the phase structure are under a slight change of 
the mass input. The second one is to get a quantitative measure for the distance 
(in mass parameter space) between the realistic and the critical meson masses, cf. 
section 5 below. As we focus on the aspect of the mass sensitivity, we have to find a 
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prescription how to tune the masses in the high- dimensional meson mass parameter 
space. On the quark level the mass parameter space is only two-dimensional, the two 
parameters are m u ^ and m s . Not only the parametrization, but also the tuning in 
the space of meson masses is by far not unique. The idea now is to parametrize the 
(pseudo) scalar meson masses by two parameters like the quark masses, the external 
fields e and e 8 . 

A relation between m u ^m s and £o,£s is obtained by identifying terms of the 
Lagrangian on the mesonic and on the quark level which transform identically under 
SU(3) x SU(3). We have 

(— e a — e 8 a 8 ) on the mesonic level, (4) 
+ {m u uu + m d dd + m s ss) on the quark level. (5) 

Thus we find 

— e = a{2m + m s ) 

-e 8 = (3{fn-m s ). (6) 

Here a and (3 are constants. They can be fixed from the known values of Eq, Eg, m Ui d 
and m s under realistic conditions. Realistic meson masses are obtained for Eq = 
0.0265 [GeV]V 8 = -0.0345 [GeV]' ', see below. The values for "realistic" current 
quark masses are taken from [HJ, rh = (m u + m^)/2 = 11.25 ± 1.45 [MeV], m s = 
205 ± 50 [MeV]. For a and (5 we then obtain 

a = -0.1164 [GeV] 2 , (3 = -0.1780[GeV] 2 . (7) 

Thus a variation in (eq,Es) can be mapped onto a variation of (rh, m s ). 

Next we have to find a mapping between (e ,e s ) and the (pseudo) scalar meson 
masses. As it is possible to explain the variety of (pseudo) scalar meson masses on 
the basis of two quark masses m u ^ and m s , it should be similarly possible to reach 
any point in an (m^, uik,...) -diagram by a variation of Eq and eg. Thus we start with 
the parametrization of the chiral limit /Xq, /1, f 2 , g, £0 = 0, £% = 0, keep the couplings 
A*o> A> hi 9 fixed and switch on e ^ 0,e 8 7^ 0. The S'f/(3)-symmetric case with 
finite, but degenerate pseudoscalar meson masses is obtained for Eo 7^ 0,^8 = 0. 
The mass point with realistic meson masses which come close to their experimental 
values, is obtained for e = 0.0265 [GeV] 3 ,e 8 = —0.0345 [GeV] 3 . (In principle, there 
is no obstacle to further optimize the values of Eq,Es such that they do reproduce 
the experimental mass values). For the chiral values /1q, fi, g and a certain choice 
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for Eq and e 8 we fi rs t determine (o~o) and (cr 8 ), the condensates at zero temperature, 
as zeroes of (§),(§) in cr and cr 8 

e + fa - -j=(2a 2 - a 2 ) + 2 -^f2°l - 4(A + |K - + f 2 )a a 2 8 = 
£s + fa* + J \gWl + vWs) + 2y/2 f 2 a a 2 8 - 4(/j + ^)ct 8 3 - 4(/i + / 2 )afc 8 = 0. 



Eqs. are the equations of motion for constant background fields ao,cr 8 . The 

pseudoscalar meson masses are then given as 



m l 



ni 



'/mi 



m 



r?08 
2 



m: 



-t4 + (4/x + ^ 2 )a 2 + (4/i + lf 2 )a 2 8 + ^/ 2 a t7 8 + -^U - 2</~ 



v/3 



14 



-/x 2 + (4/i + ~/ 2 )<7 2 + + - ^V2f 2 a (7 8 + ^=ga + \ - 



-/io + (4/i + + (4/x + ^/ 2 )a 8 2 - ^=ga 

-Vl + (4/i + \hK + (4/j + 2/ 2 )a 8 2 - ^V2f 2 a a 8 



08 



8 S 2 r- 2 

- gV2/ 2 (j 8 + -j=g 



\( m l m + m is - \/Koo - m L) 2 + 4 Ko8) 2 

O (""4, + ^SS + \/KoO - m L) 2 + 4 ( m 



2 )2) 
'708 '' ' 



V 7 ^ 



0^8 



v/3 



9<T0 + 2\ -g^S 



(10) 



The scalar meson masses follow from 



m 



n? 



coo 



m 



in 



C08 
2 



"7 



-/i 2 + (4/! + 4/ 2 )a 2 + (4/i + 2/ 2 )a 2 + ^=f 2 a a 8 - -=ga + 2\ -ga 8 



V2 J 



-(m 2 oo +m 2 88 + ^(m 2 
1 



CT00 -< 8 ) 2 + 4« 8 ) 2 ) 



2 , m a 00 + m *88 ~ V K» - m2 8 8 ) 2 + 4(^ 2 08 ) 2 ^ 



;g<?8 



-/i 2 + (4/ x + 4/ 2 )a 2 + (4ft + 2/ 2 )a 2 - 4^ Ws - ^9*0 ~ y o 
-/i 2 + (12/x + 4/ 2 )a 2 + 4(/i + / 2 y 8 + i|a 

-/i 2 + 4(/i + / 2 K 2 + 6(2/! + / 2 )cr 8 2 - 4 ■ V2/ 2 oo<t 8 - -j=v - 2^|po- 8 

4 2 
o" 8 [8(/i + / 2 )o- - ^7^/2^8 - -j=9] 



(11) 
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where a , <r 8 should be replaced by (cr ), (cr 8 ), respectively, everywhere in Eqs. (10, 11) 
for physical masses. The resulting masses together with the input parameters are 
listed in Table 1 of section 6. 



To facilitate a comparison with other work on the chiral phase transition, it 
remains to translate the results for the meson condensates to the quark level. If we 
identify terms in the Lagrangians on the quark level and on the mesonic level (cf. 
Eqs. (g) and (||)), we find 

\{2rh + m s ){2qq + ss) = -£ o"o (12) 
2 

-(m - m s )(qq - ss) = -e 8 a 8 (13) 
with qq=\ [uu + dd) , leading to 

^ = 2^^ ( " o) + 2(m"-m s ) ((T8) 

(ss) = ^^-M + ^-fo). (14) 
2m + m s m — m s 

The coefficients are just proportional to a and ±/3, as a comparison with Eq. (^|) 
shows, a and j3 have been determined in Eq. (0). Later we take the relations ( fU]) 
as temperature independent and substitute (o"o)(T), {a%){T) for the corresponding 
condensates at zero temperature. This assumption is consistent with our approach. 
We also treat the couplings /ig, /i, f2,g of the Lagrangian as temperature indepen- 
dent. Therefore the symmetry of L remains unchanged under an increase of T. The 
symmetry was on the basis of the identification which has led to Eqs. (|I4D. (The 
assumption of temperature independent couplings may not be justified in the vicin- 
ity of the transition region.) In the next section we will outline how to calculate 
(a >(T)> 8 >(T). 



3 Large A^-expansion 

In an earlier calculation the linear SU(3) x SU(3) sigma-model has been considered 
in a mean-field approximation ||20|| . Recently, Gavin, Goksch and Pisarski have 
tried to localize the first order transition boundary in an {m u ^, m s )-mass diagram in 
a mean-field calculation. The famous renormalization group analysis of Pisarski and 
Wilczek applied to the linear sigma-model in the chiral limit. Frei and Patkos ]2T 



were the first to apply a saddle-point approximation to the partition function of the 
sigma-model. Their investigations were also restricted to the chiral limit. Meyer- 
Ortmanns, Patkos and Pirner 0] have extended the approach of Frei and Patkos to 
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finite meson masses. In [|TTJ only the zero-Matsubara frequencies were kept. When 



the imaginary time dependence of the fields or, equivalently, the n 7^ O-Matsubara 
frequencies are dropped, it results in a dimensional reduction of the four- dimensional 
theory to an effective three-dimensional theory. In general, such a reduction can be 
justified in the high-T-limit or for an anticipated second order phase transition. In 
both cases the ratio of is negligible (/? is the inverse temperature and £ denotes 
the largest correlation length of the system). The n 7^ O-Matsubara modes were 



included in ||18|| . In this paper, we follow the same approach as in [18], but extend 



the work to study the aspect of the mass sensitivity of the chiral transition. Further 



differences in details of |L8[ and the present paper will be mentioned below. In the 
following, we summarize the essence of the large-iVy-approach to make the paper 
self-contained. 

The temperature-dependent order parameters are the meson condensates (<To)(T) 
and (ag)(T). They are related to the light and strange quark condensates according 
to Eqs. fllip . The values of (o"o)(T), (as)(T) are determined as the minima of an 
effective potential U e ff(ao, cr 8 ). The effective potential is calculated as a constrained 
free energy density, i.e. the free energy density under the constraint that the average 
fields -j^y Jq dr J (PxaQfi(x,T) take prescribed values cr and cr 8 , respectively, while 
the same averages over 07, £ = 1, 7, and ire, £ = 0, 8 should vanish. The phys- 
ical meaning is obvious. If one chooses Oo^s 7^ in the high-temperature chiral 
symmetric phase, the corresponding free energy is certainly not minimal for such a 
choice. We consider spontaneous symmetry breaking in two directions. Accordingly, 
we introduce two background fields <r , a 8 

cr 8 = o-g + o-g, (15) 

where a' Q , a' 8 denote the fluctuations around the background. Otherwise, oi = a' e 
for £ = 1,...,7 and 7T£ = 7r' e for £ = 0,...,8. As a common notation for cr^rr'e we 
use M' e = <j\ + m' t . The actual minima of U e ff(ao,cr 8 ) (with a = o"o,°"8 = ^s) are 
denoted as (a ) and (cr 8 ). 

The Lagrangian is expanded in powers of M[. The linear term in M[ vanishes due 
to the (^-constraint in the constrained free energy density. The quadratic term defines 
the masses ttiq of the meson multiplets ir, K, 77, 7/, a n , <tk, oy , Q = 1, 8 labels 
the multiplets. The explicit formulas were given in Eqs. (|T0|) , ([TT|) . The two quartic 
terms are quadratized by introducing an auxiliary matrix field XX X ) according to 
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2j 

exp{-/3[/ 1 (TrM'M /+ ) 2 + / 2 Tr(M' M' + ) 2 ]} = 

(16) 

rc+ioo 

= const ■ / DS(x) exp{Tr£ 2 + 2eTr(EM'M /+ ) + 2aTr(M'M')TrE}, 

•/ c—ioo 



' C— lOO 

where M'(x) is an iV x iV-matrix field and 



e 2 = Pf* 



2ea + 3a z = Pf x . 

Note that Eqs. ([U]) are a sophisticated version of the simpler case, where is a 
scalar field 

const ■ exp{-a(0 2 (x)) 2 } =/ PS(i)e 22w - s(l) ^ (x)2 ^. (17) 

Formula (|PT|) can be easily generalized to the case, where the l.h.s. includes also a 
cubic term <f) 3 (x), but we are not aware of an analogous transformation that leads 
to a tractable expression if the cubic term occurs in the form of a determinant. 
This is the reason why we drop the cubic term in M' e in the following procedure. 
The (3 x 3)-matrix field E(x) is replaced by an S , f/(3)-symmetric diagonal matrix 
£ = diag(sa<i, sad, sad). Thus the matrix of auxiliary fields is reduced to a single 
field variable sad(x). The quadratization of the quartic term of the Lagrangian 
leads to 



£(4) = / 1 (TrM'M /+ ) 2 + f 2 Tr(M'M 



>+\2 



- £ (4)/ = 3 - {^-sad 2 + faad) +\{s + /. 2 )Tr(M'M'+). (18) 

In the saddle-point approximation the path integral / PE(x) is dropped, the 
auxiliary field sad(x) is replaced by sad*, which maximizes the integrand. For an 
0(iV)-model it is well known that this approximation corresponds to the leading 
order in a 1/iV-expansion [^|. In the special case of f 2 = = g, the SU(3) x SU(3) 
sigma-model becomes an 0(18)-model, which is invariant under 0(18)-rotations. 
We have N = 18 = 2iV 2 , where Nf denotes the number of quark flavours, while 
N labels the number of mesonic modes. Terms of 0(1/N) are dropped, as long 
as fluctuations in the auxiliary field are neglected. Therefore we call our scheme a 
large — Nf — approximation. 

The resulting one-loop contribution to the free energy density is given as 
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where g(ir) = 3,g(K) = i,g(r}) = l,g{r]') = l,g(a n ) = 3,g(cr K ) = 4,,g(a v ) = 
l,g(a v ') = 1 are the multiplicity factors of the multiplets, uo 2 n = (27m/T) 2 and 

X 2 Q = sad + fx 2 + m 2 Q . (20) 

Thus the one-loop contribution to the free energy takes a form, which is familiar from 
a free field theory. The only remnant of the interaction is hidden in the effective mass 
square Xq via the auxiliary field variable sad. After the sum over the Matsubara 
frequencies is performed, the full expression for the effective potential is given as 

U ef f(a , a 8 , sad) = U c i ass (a , er 8 ) + U sadd i e (sad) + U (a , cr 8 , sad) + U th (a , a 8 , sad), 

(21) 

where a = a , a$ = a 8 and 

U dass (a , a 8 ) = + vl) + ^( 2 ^o - ^4 ~ ^l) (22) 

2V / 2 X 3 i / t , /2\ 4 , ft , fas & . n ft , f 



■hw\ + {h + + (/i + -KM + 2(/i + f2K<4 ~ Wo - e 8 a 8 



3 ° w 3 ' u w 2 

is the classical part of the potential, which is independent of sad. 

U saddle {sad) = - 8( ^ 3 + {~^~ + l4 sad^j (23) 

results from the transformation ([16]), cf. Eq. (fL8|) . The one-loop part consists of the 
zero point energy Uq and the thermal part Uth, given by 



1 /" A d 6 K 

U (a , a 8 , sad) = g(Q) J j^y^Qi (24) 



t/p is divergent if the three-momentum cut-off A is sent to infinity, while U t h is 
convergent for A —>■ oo and vanishes for T = 

^4g s(o) /S in(i - e ^ o) - (25) 

Our results have been obtained for the potential U e ff = U e ff — Uq, where the diver- 
gent zero-point energy has just been dropped. It is often argued that the omission 
of the zero-point energy is justified if one is finally interested in thermodynamic 
quantities, which are derived as derivatives of InZ w.r.t. (3. Usually, Uq does not 
depend on T, and the splitting of the one- loop part of U e ff according to Uq and 
U t h is a splitting in a T-independent and T-dependent part. Thus the contribution 
of Uq should have no effect on the thermodynamics. Strictly speaking we cannot 
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use this argument, because U of Eq. (|24|) has an implicit temperature dependence 
hidden in Xq via sad(T), the temperature dependent saddle point variable, which 
is finally chosen such that it maximizes U e ff for sad = sad*. This was one of the 
reasons why the [/o-term was kept in [|I8| . A renormalization prescription was im- 
posed such that the strong (quartic) cut-off dependence of U was weakened to a 
l/A 2 -dependence. Further differences to |18[ are due to corrections of two errors 



in flBfl . After we had removed the programming error in U c i ass of [F^], we found 
two branches in the free energy density / (or the pressure). The branches in / of 
the high and low-temperature phases did not cross at some temperature T = T c , 
although the free energy should behave as a smooth function in T by general ther- 
modynamic arguments. The reason was that the two minima of the potential which 
exist below T = 200 [MeV] were erroneously associated to the true minima in the 
high- and low-temperature phases. It was not realized that for all temperatures up 
to T < 200 [MeV] one minimum is only a local one, while the other one stays the 
absolute minimum for all temperatures up to T > 200 [MeV], cf. also the remarks 
in section 6. 

Thermodynamic quantities. Thermodynamic observables are derived from the 
free energy density / = limy^oo(— -Jy In Z), where Z is approximated as 

Z«<7 >, fa)) = e -^^/((-o>(T),<<x 8 >(T) ;s ^(T))_ (26) 

We have explicitly indicated that U e ff should be taken at the saddle point value 
sad*(T), which extremizes U e ff- The expression for the energy density e, entropy 
density s, and the pressure p are derived in the standard way. We have 

V = -Ueff, (27) 



I "* poo 

L giQ)i ^Jo da 



813 tt»^ J 2ir*h \ P V^Q 



a 



K/T, y Q = (yW + ^ + m^ao),^))) /T. (29) 



Evaluation of U e ff. Next let us discuss the evaluation of U e ff m some more 
detail. The expression for U e ff is familiar from the expression for a free field theory. 
In a high-temperature expansion it reads 

Ueff = U dass + U saddle -|^E 9(Q) - ^Vq + 0(y 3 Q ) | (30) 
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with y Q = 



At high temperature the SU(3) x SU{3) linear sigma-model certainly fails to 
describe the quark-gluon plasma phase. The critical temperature falls neither in a 
high nor in a low-temperature region, the expansion parameter is of 0(1) or larger. 
Nevertheless, we discuss the high-temperature expansion, which has been frequently 
applied in calculations of a thermodynamic potential. It has the advantage that 
imaginary parts of U e ff are absent to leading order in Xq/T, cf. e.g. p3] . We 



have performed a high-temperature expansion for two sets of mass parameters, the 
chiral limit and the realistic mass point, see sections 4 and 7. This way we got 
some qualitative insight in the phase structure before we could tackle the problems 
in a fully numerical evaluation of U e ff- Analytic expressions for a low-temperature 
expansion of U e ff are known as well |24|. We have used them only in intermediate 



steps to check the numerics. Thus we turn now to the fully numerical evaluation of 

Ueff- 

At a first glance the numerical evaluation of U e ff looks rather straightforward. 
For each pair of (o"o, erg) we have to find sad* that maximizes U e ff. The condensates 
(o" ), (erg) are then determined as minima of U e ff(a , cr 8 ; sad*(a , <7 8 )). It is well 
known [23] that the arguments (K 2 + Xq) of the logarithm in U e ff can become 



negative and lead to imaginary parts in U e ff, which we have mentioned above. The 
original hope was that the positive contribution of the auxiliary field sad to the 
masses tuq helps in avoiding imaginary parts of U e ff. Actually no imaginary parts 
have been found in |?TJ. As can be seen from Fig. 2 the contribution of sad* to 



Xq increases with T, it is positive for T > 116 [MeV], but in our case the positive 
contribution is not sufficient to avoid negative Xq completely, i.e. for all <r , <t 8 and 
T. It has turned out that the actual minima (ctq), (<jg) lie always in the "allowed" 
region of real valued U e ff or at least at the boundary of this region, but on the way of 
searching the maximum in sad* and the minima in (cr ), (cig) the routines encounter 
negative mass squares indispensably. Therefore U e ff is analytically continued in the 
following way. The integral in U th of U e ff is of the type 

/(*)=/ dx{x 2 lull - exp\-{x 2 + z 2 ) 1/2 }} = -z 2 y —K 2 {nz), (31) 
Jo ' ,71 " 

where K 2 is a modified Bessel function. The analytic continuation of K m (nz) to 
complex values of nz is given as 



K m (nz) = i 



m+l 



(0 [J m (inz) +iY m (inz)), (32) 



where J m are Bessel functions of the first kind and Y m are Weber functions 25 
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Figure 2: Saddle point sad* as a function of temperature T 



This leads to the following form of U e ff 

T 4 8 ( —X 2 00 1 / n 

Ueff = U dass + U saMle + — ^9(Q)^^^^K 2 (^-X Q )\ (33) 



for > 0. For X 2 Q < and X Q = ±iJ|Xg| 

#e// = U dass + U sadd i e + Re U th + i\vsiU th 

with 

rpi OO f I ~V^2 CO 1 / s ^ 

^ - |^«){±^£>(^r)}- 



(34) 



In the actual calculation we have cut the infinite sums in Eq. after 25 

terms. 
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The condensate values (<r )(T), (<r 8 )(T) are determined as minima of U e ff for Xq > 
and of Re U e ff for Xq < 0. 

Error analysis. There are mainly two sources for errors. The first one is the 
determination of the saddle point from the condition dU e ff/dsad = 0. In searching 
sad* the routines encounter negative values of Xq in intermediate steps and find sad* 
from derivatives of Weber functions Y2(^\Xq\ 1 ^ 2 ), see Eq. (|34|). The Weber func- 
tions are taken from representations in [^6j which cover the range from small to large 
arguments. For arguments ~\X%\ 1/2 > 4.8(9.7), Y 2 (^\X 2 \ 1 / 2 ) < 6.0 • l(r 3 (2.9 • 1(T 5 ). 
Due to the large variation of the effective masses the arguments become easily > 100. 
Thus for small temperatures, large effective masses and/or large values of n the con- 
tributions to the sum over n are so small that they reach the order of the numerical 
accuracy. The location of the saddle point becomes inaccurate under these condi- 
tions. Similarly, for small arguments (i.e. high temperatures, small values of Xq) 
the approximations of the Weber functions in the vicinity of their singularity at 
vanishing arguments become less reliable p6| |. This explains why the error in sad* 



increases with temperature, if simultaneously jXgl 1 / 2 becomes smaller. 
The error Asad* in the maximum sad* of U e ff is estimated as the difference between 
the zeroes of dU e ff/dsad = if sad* is extrapolated from the right or approached 
from the left in the following way. On the left hand side of the maximum, Re U e ff 
is represented by the series of Weber functions, it is oscillating and dU e ff/dsad > 0. 
The zero of the oscillating derivative is actually used in the program as value for 
sad*. The oscillations extend to some region (of temperature dependent size) on 
the right hand side of dU e ff/dsad = 0, where the derivative is negative. For further 
increasing values of sad the effective potential becomes real and behaves approxi- 
mately quadratically in sad, such that dU e ff/dsad is linear in sad. We extrapolate 
the linear part of the derivative to zero and use the extrapolated zero as second 
value for sad*. The difference between both values for sad* gives an estimate for 
the ambiguity in finding the maximum of U e ff from the parabolic or the oscillating 
behaviour. It turns out that for low temperatures Asad* ~ 0, while it increases 
with temperature to Asad* ~ ±0.006 [GeV] 2 in the transition region in case of the 
chiral limit, cf. Fig. 2. 

The uncertainty in the determination in sad* leads to errors in the meson conden- 
sates, the effective masses Xq, and the thermodynamic quantities e,p, and s. The 
errors A{ao^) which are induced in (cr ), (cs) due to the uncertainty in the precise 
value of sad* are estimated from 



- d 2 U ef f/dsad* da 



d 2 U eff /d 2 a 2 



■ Asad*, (35) 

(cr ),{a s ),sad* 
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since (a 0)8 ) (sad*) are implicitly determined via dU e ff/da 0t8 = 0. In case of realistic 
masses this leads to A(cr , 8 ) < 10~ 7 below T ~ 50 MeV, < 10~ 4 below T ~ 125 MeV, 
< 10~ 3 for T < 165 MeV and of the order of 10~ 3 [MeV] in the transition region. 
After the crossover (T > 197 MeV) it stays ~ O(10~ 4 ). 

An independent second source for errors in the condensates (cr ), (erg) is the flat 
shape of the effective potential in the transition region, in particular for external 
field strengths which only admit a rather weak first order phase transition. We 
estimate A((j )(— A(crg)) from plots of £7 e //(a"o) in the 5'L r (3)-symmetric case. We 
find A(ct ) - ±1-3 [MeV] in the chiral limit and ~ ±4.3 [MeV] for e = 2 • 10~ 4 or 
2.5 • 10~ 4 [GeV] 3 , cf. Fig. 6 of section 5. 

In the S{7(3)-symmetric case (as) should identically vanish for all temperatures. 
The eight Goldstone bosons are identical in mass. If the S"?7(3)-symmetric case 
is calculated as a special case (eo ^ = e%) of our general framework (e ^ ^ 
£g), (as) = should come out automatically. In fact it vanishes within the numerical 
accuracy for temperatures T < 50 [MeV] for eo = 2.5 -10~ 4 [GeV] 3 . In the transition 
region (a$) strongly fluctuates around zero. The maximal fluctuation varies between 
0.2 and 3.0 [MeV] depending on e and the strength of the first order transition. 
The values are of the same order of magnitude as A(ao,s) due to the flat shape of 
the effective potential. 

When the meson condensates are converted to the quark condensates, additional 
errors enter due to the current quark masses. We use Am u ^ = ±1.45 [MeV] |i~9f . 
Finally, the errors in the entropy and energy densities s and e can be traced back 
to the effective masses Xq((ao), (a s ), sad*). The resulting errors in s/T 3 in the 
crossover region are indicated in Fig. 10 of section 6, where s/T 3 is shown for realistic 
meson masses. 

From very different approximation schemes it is known that the region of a phase 
transition is particularly difficult to handle. We consider our difficulties in localizing 
the saddle point and calculating the condensates in the transition/crossover region 
just as a special manifestation of that. 

4 The chiral limit 

From the renormalization group analysis of Pisarski and Wilczek we expect a 
first order chiral transition in the chiral limit. As mass input we choose = m 2 K = 
m 2 = for the pseudoscalar Goldstone bosons, = 850 [MeV], m a , = 800 [MeV], 
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f n = 94 [MeV] . The mass of the sigma meson should be treated as a parameter due 
to the uncertainty in the experimental identification. For the final parametrization 
we choose m ari = 600 [MeV]. 

The tree level parametrization of the Lagrangian ([I]) in the limit of vanishing e 
and e 8 is then given by 

l4 = 5.96 • 10" 2 [GeV] 2 
fi = 4.17 
h = 4.48 

g = -1.81[GeV]. (36) 
The high-temperature expansion gives a 1st order transition at T c = 92 ± 1 [MeV], 



Figure 3: Normalized light quark condensate (qq)r/ (qq)o as a function of temper- 
ature T in the high-temperature expansion (solid curve) and the numerical evalu- 
ation (dashed curve) of the large iVy-expansion. The condensate drops to zero at 
T c ~ 92 [MeV] and 177 [MeV], respectively. 

while the fully numerical evaluation leads to T c = 177±1 [MeV], cf. Fig. 3. This is in 
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agreement with the general observation that the high-T-expansion gives qualitatively 
correct results when it is extrapolated beyond its validity range, quantitatively it 
fails in precise predictions. The T c -value of the numerical calculation supports the 
estimate for the iV/- dependence of T c , which has been derived by Cleymans, Kocic 
and Scadron |27j from a pion gas model without interactions as 



T c ~ 2 • U ■ ^N f /(Nf - 1) ~ 200 [MeV]. (37) 

The effective potential U e ff is plagued with imaginary parts for all temperatures 
we have investigated, i.e. between and 250 [MeV]. When we evaluate Re U e ff 
according to Eq. fl3"4[) , we observe the same effect for higher values of the sigma- 



mass as it has been noticed by Goldberg p0 |. Goldberg has calculated the effective 



potential of the SU(3) x SU(3) sigma-model in a mean-field approximation in the 
chiral limit. The oscillations in the potential became stronger for larger values of the 
cr-meson mass, which was used as input. In our case Re U e ff(a ) strongly fluctuates 
around the expected parabola if m a > 1 [GeV]. This explains our final parameter 



choice of m ari = 600 [MeV]. The approximation of Re U e ff according to Eq. fl3~4]) by 
a finite series of Weber functions (n < 25) loses its validity if the argument of Yi 
becomes too large due to m|. 

For the barrier height between the coexisting minima of the potential we find 0.14 • 
10 -3 [GeV/fm 3 ]. The barrier height may be regarded as a measure for the strength 
of the transition. It determines the tunneling rate between coexisting phases in the 
transition region. The barrier height is clearly smaller than the value of Frei and 
Patkos [^TJ, who found 0.36 [GeV/fm 3 ] in the same model, but without inclusion of 



the n 7^ 0-Matsubara frequencies. The higher barrier goes along with a large value 
for the interface tension a between the coexisting chirally broken and symmetric 
phases at T c , a has been estimated as [(40 — 50)MeV] 3 More recent lattice 

results indicate that such a large value is not likely for QCD 



Of particular interest in the chiral limit is the temperature dependence of the 
effective masses Xq = sad + /ip + rrin, Q — 1, ...,8. In the chiral limit the masses 
are free of contributions of the external fields. Thus the Goldstone theorem should 
be obeyed by the pseudoscalar meson octet. A short glance at Fig. 4 shows that 
the effective masses Xq = X v = X^ = X n do not stay massless for T > 0. They 
get an increasingly positive contribution with increasing T due to the saddle point 
contribution sad*(T). If the physical meaning of Xq coincides to leading order in 
1/Nf with a dynamical mass (in particular if the dynamical and screening masses 
coincide to leading order in 1/Nf and both are equal to Xq), the increase of Xq(T) 
with T would violate Goldstone's theorem and should be considered as an artifact 
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Figure 4: Effective masses Xq , Q — 1, 8 as a function of T in the chiral limit. 
In units of [MeV] we have at T = X x = X K = X v = 0,X V , = 850, X a , = 600, 
= X aK =X av = 800. The masses are degenerate at T c ~ 177 ± 3 [MeV]. 

of our approximation. We still have to clarify this point. 



5 The critical transition line 



The critical transition line in an (m n , rriK, ■ ■ .)-diagram consists of pseudoscalar me- 
son masses for which the first order chiral transition becomes of second order and 
turns into a crossover phenomenon for meson masses exceeding the critical values. 
We will determine three such critical points. The critical points are characterized 
by their external field strengths eo,es- The first point corresponds to an SU(3)- 
symmetric case, where e$ = 0, m u = = m s ^ and m n = itlk = vn^ ^ 0. 
Since < a 8 >= for all temperatures, the numerics considerably simplifies com- 
pared to the general case with eq ^ ^ eg. For the second critical point we choose 
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6q = —0.77 eg. This ratio is identical to Eq/es for the mass point with realistic 
meson masses, where m s /rh = 18.2, cf. section 2. The third point is characterized 
by e = (2a/ P) ■ £§• It is chosen such that m s = 0, m ud ^ 0. Before we present 
our results for the critical field strengths in the large-iVy-expansion, we calculate 
£q"*, E^ lt in a mean-field approximation. 



5.1 Critical meson masses in a mean- field approximation 



Recently, Gavin, Goksch and Pisarski 17 have calculated a set of critical quark 



masses in the linear SU(3) x SU(3) sigma-model. The calculation has been per- 
formed in a mean-field approximation. Since we use a different tree-level parametriza- 
tion of the sigma-model, we have reproduced the mean-field calculation for our pa- 
rameter choice and summarize the main steps. 

The SU (3) -symmetric case. Let us first consider the S*?7(3)-symmetric case with 
£8 = 0. In mean- field we have to deal with the classical part of the potential U c i ass , 
which follows from the Lagrangian flU) for a constant background field do = a (ct 8 = 
0). We have 

12 A 
UdassM = --/io^o + ^"y= #°o + (fi + y Vo ~ £ o°o- (38) 

This form is familiar from Landau's free energy functional in terms of an order 
parameter field. One of its applications is a description of the phase structure for 
a liquid/gas transition. For e = the system has a first order transition from the 
liquid to the gas phase. The transition stays first order, until the external field (the 
pressure in case of a liquid/gas system) is increased to a critical value e^, where 
it becomes second order. For values of £q > E^ rit , the transition is washed out and 
turns into a crossover phenomenon. 

In mean-field the effect of a finite (and strictly speaking high) temperature is 
reduced to a renormalization of the mass parameter term, i.e. of the coefficient of 
the quadratic term in the Lagrangian. Thus the finite temperature effects can be 
mimiced by tuning /Xq while keeping the other couplings in the Lagrangian (fx, f 2 , g) 
fixed. The potential can be Taylor expanded around its true minimum a™ in (which is 
different from zero in the symmetry broken phase); in particular it can be expanded 
around the "critical" a™ m = erg"* for the critical field strength e c q %1 . In the "critical" 
case U c i ass starts with a term quartic in (ctq — e^™*). From the vanishing of the first 
three coefficients we obtain the critical parameters as follows. 
Since the first order transition just disappears for Eq" 1 , the third derivative of U, 



class 
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w.r.t. ctq should vanish at o~ — Cq leading to 



a cr it = 9 = 3 A _ 10 -2r GeV ]. (39) 

The second derivative d 2 U c i ass / 'da 2 at Cg™* is the coefficient of the quadratic fluctu- 
ations (do — ctq"*) 2 around ctq™*, thus it has the meaning of m 2 r This is the critical 
mass, which goes to zero when the second order transition is approached. It is eas- 
ily checked that the other meson masses remain finite for the same set of critical 
parameters. The vanishing of m 2 t at criticality or d 2 U c i ass /da 2 at erg"* implies for 
the critical value of /i 2 , 

£ «* = Q ; m = -6.44 ■ 10" 2 [GeV] 2 . (40) 



Finally the extremum condition for a^ %t determines e^ lt as 

i g 3 

'27-6-VW1 + /2/3) 2 



4 rit = 1 rit —. = 6.6 ■ 1(T 4 [GeV] 3 . (41) 



The value for e c ^ %t = 6.6 • 1CT 4 [GeV] 3 (£§"* = 0) leads to a pseudoscalar octet 
mass of m n = itlk = rn v = 146 [MeV] and to a current quark mass of m u = 
m d — m s — 1-9 [MeV]. One should keep in mind that the tree level parametriza- 
tion of the SU(3) x SU(3) sigma-model is arbitrary to some extent. If we would 
choose g = -1.39 [GeV], ft = 5.3, f 2 = 0.93 with m % , = 600 [MeV] (the val- 
ues which have been used in JTT| for the tree level parametrization), e^" lt turns out 



as 3 ■ 10~ 4 [GeV] 3 or m uAs ~ 0.9 [MeV] leading to (m ps ) = 115 [MeV]. For the 
same parameter choice, but f\ = 2.35 and m ar] = 950 [MeV],£Q rii comes out as 
1.3 ■ 10~ 3 [GeV] . The same tendency has been observed in |T7|]. An increase of 
the input mass m a shifts the critical field strength to larger values (reducing the 
discrepancy to the lattice result). The parameter choice which has been used in 



leads to £q"* = 4- 10~ 4 [GeV] . Gavin, Goksch and Pisarski [17] obtain for the corre- 
sponding critical field strength in the 5'?7(3)-symmetric case h c ^ %t = 1.6 ■ 10 -4 [GeV] 3 . 
Due to their different parametrization it is not identical with ours, but of the same 
order of magnitude. Thus we estimate the error in the mean-field calculation as 
Ae° rit = ±5 ■ 10~ 4 [GeV] . The induced error in the pseudoscalar meson masses 
comes out as large as ~ 127 [MeV], in the current light quark mass it is < 1.4 [MeV]. 



The general case of e ^ 0, 0. For non-vanishing e and e 8 U c [ ass is given 

by Eq. (p2[). Again we have to determine the critical parameter /zjj to simulate a 
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finite temperature, the critical minima values ctq 7 " 2 *, a 8 rit of U c i ass and the critical 



field strengths e 



crit 



-crit 



We use the following conditions 



dUd ass (ao,a 8 )/da 
dU dass (a ,cr 8 )/das 



m„ 



e^/e c 8 rit + 0.77 



crit 



crit 



crit 



crit 



d 3 U c i ass (cr , <7 8 )/dr 3 \ crit 











The five conditions are postulated at criticality (abbreviated as 



(42) 
(43) 
(44) 
(45) 
(46) 

i.e. for the 



crit 



r crit 



e C Q %t : , s 8 %t , /1q cnt . The first three equations 



set of critical parameters (Tq, a 
are obvious generalizations of the SU (3)-symmetric case. Eq. fl4"o] ) is just one pos- 
sible choice saying that we keep the mass splitting of realistic (pseudo) scalar me- 
son masses fixed in the tuning to the critical transition line. Eq. (^) generalizes 
d 3 U c i ass (ao) / 0(Jq \ crit = of the one- dimensional case. More precisely, it postulates 
that the directional derivative in the radial direction r in (<7o, crs)-space (i.e. per- 
pendicular to the direction of m 2 = 0) should vanish to exclude the occurrence of 
a first order transition. Eqs. (|42| - ^3|) lead to 

2 



.crit 



.crit 



9 



(2a - a 8 ) - -V2f 2 a 3 + 4(/x + f 2 /3)a + 4(/j + f 2 )a a 2 8 



y/2a cr 8 ) — 2y/2f 2 aoul + 



+4(/ 1 + / 2 /2)a 8 3 + 4(/ 1 + / 2 )a V 8 , 



(47) 



crit 



) a 8 



a 



crit 



where finally Hq = /Zq cnt , cr = a, 
If we introduce the auxiliary quantities ro,r 8 according to 

r = (12/x + 4f 2 )a 2 +4(f 1 + f 2 )a 2 8 +Aga /V3 

implies 



4(/x + h)°\ + 6(2/! + f 2 )al - 4V2f 2 a a 8 - ^=a - 2^ga 8 (48) 



and use m^ og of Eq. ([□]) , Eq. 

1 
2 

where in the end o~o = er. 



2 crit 



r + r 8 ~ J ( r o - r s) 2 + 4(m2 Qg ) 2 



crit „ 
; a 8 



of"* in Eqs. (g8| 



So far we have e^g* 



£ o!8*(/i> /a> 9, o"o, crs) and /xg 



2 crit 



,,2 crit 

Mo 



(49) 



{fi,f2,9,<ro,<r%) 



Finally the zeroes of Eqs. ( |4ll| - f4"6|) determine cr^™ 4 , a 8 nt . These equations are solved 
numerically. We find 



a™ = 3.6- 10^[GeV], a c s 



-5.5- 10" d [GeV]. 



(50) 
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Using these values, Eqs. (g9j), (f47|) and (|48|) give 

^crit = _ 6 .5 . i - 2 [GeV] 2 , (51) 
e% lt = 7.7- l(T 4 [GeV] 3 , ef u = -1.0 • l(T 3 [GeV] 3 (52) 

or m s = 6 [MeV], m = 0.3 [MeV], while the average pseudoscalar meson mass {m ps ) 
= 126.4 [MeV]. The values for are compatible with the results of Gavin, Goksch 
and Pisarski |T7J, who find for the critical field strengths h^ %t = (62 [MeV]) 3 , h^ %t = 
(60.4 [MeV]) 3 , if one keeps in mind the different tree level parametrization. For 
example m s /m Uj d = 32 in |T7[, while m s /m u ^ = 18.2 in our case. 

The case ofm s = or eo/e$ = 2a/ (3. If we replace 0.77 in Eq. (|45l) by — 2a/ (3 = 
— 1.31 corresponding to m s = 0, m u ^ ^ 0, we obtain from Eqs. ((421) -(flB]) 



al nt = 3.4- 10- 2 [GeV], a% u = 2.9 ■ 10~ 3 [GeV] (53) 
^crit = _ 6 . 4 . i - 2 [GeV] 2 (54) 
eg** = 6.7- 10" 4 [GeV] 3 , 1 = 5.15 ■ 10" 4 [GeV] 3 (55) 

or rh = 2.9 [MeV] and (m ps ) = 137.2 [MeV]. Note that the condition m s = implies 
the same sign for e^" 11 and ef lt due to the same sign for the constants a and (3. 
Hence the "critical" condensates cr^g* come out with equal sign in contrast to the 
realistic mass case. 

In the following we will compare the mean- field values for e^f for Eqs. (41), (17). 
((5"5|) with the large- iVy-results. 



5.2 Critical meson masses in the large- iVj- expansion 

The SU (3) -symmetric case. When e 8 = and Eq is slowly increased from to 
~ 2.5 • 10 -4 [GeV] 3 , we observe a weakening of the first order transition as is seen 
in Fig. 5. For the critical field strength we find 

e' rit < (3 ± 0.5) • 10- 4 [GeV] 3 . (56) 

The value for e^ lt induces a pseudoscalar meson mass of 51 [MeV] or quark masses 
of m£2f = m c ; il = 0.9 ± 0.14 [MeV], such that 

™Z!tJrn uAs ~ 0.08 ±0.01. (57) 

Within the errors the result for e^ lt is of the same order of magnitude as the 
mean-field value. The uncertainty in our result is Ae^ lt < 0.5 ■ 10~ 4 [GeV] 3 . For 
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Figure 5: The light quark condensate normalized to its value at zero temperature 
(qq)r/ (qq)o as a function of T in the S , t/(3)-symmetric case. The weakening of the 
first order transition is obvious, when e [GeV] 3 is varied between e = (x),2 • 
10~ 4 (O), 2.5 • 10" 4 (+) and 6.6 • 10" 4 (□). 

e = 2.5 • 10~ 4 [GeV] 3 the transition could be clearly identified as first order from 
the (qq)T /(qq) o(T)-curves. For e = 3 • 1CT 4 [GeV] it is a crossover. We could have 
further improved the accuracy of e§ %t by measuring data for (qq) T (T) in the inter- 
mediate 5 - ran g e - On the other hand such an improvement is anyway limited by 
the well-known fact that it is in general hard to disentangle a very weak first order 
transition from a rapid crossover phenomenon. 

The weakening of the first order transition is also revealed in the barrier height 
of the effective potential between the coexisting minima at T c . The barrier decreases 
from 1.4 • 1CT 4 [GeV/fm 3 ] in the chiral limit to 2.1 • 1CT 6 [GeV/fm 3 ] for e = 2.0 • 
10~ 4 [GeV] 3 , 6s = [GeV] 3 , the largest value for which a first order transition could 
be identified from the shape of the effective potential, cf. Fig. 6. A comparison 
between Figs. 5 and 6 shows the ambiguity in identifying a very weak 1st order 
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transition. Fig. 5 suggests a weak 1st order transition for e = 2.5 • 10~ 4 [GeV] 3 with 
T c ~ 181 [MeV], but no barrier is visible between the two coexisting condensate 
values at the same temperature and the same e - va hie in Fig. 6. Accordingly a 
precise determination of T c is hard in case of a weak first order transition such that 
we estimate for the error in T c AT C = ±2 [MeV]. Fig. 6 also admits an estimate of 
the error in finding (a 0t8 ) in the transition/crossover region if the effective potential 
is very flat, which we have mentioned in section 3. In the chiral limit we have 
A(cr ) ~ 1.3 • 10 -3 [GeV], for e = 2 • 10 -4 and 2.5 • 10" 4 the error is estimated as 
4.3 • 10~ 3 [GeV]. Later we assume A(cr ) ~ A(cr 8 ). 



Figure 6: Decrease of the barrier height of the effective potential U e ff for e = 
at T ~ 176 [MeV] (solid curve), e = 2 • 1(T 4 [GeV] 3 at T ~ 180 [MeV] 3 (dashed 
curve), £o — 2.5- 10~ 4 [GeV] 3 for T ~ 181 [MeV] (determined from the corresponding 
condensate curve of Fig. 5) (dotted curve) in the S*?7(3)-symmetric case. 

T c (m)- dependence in the SU '(3) -symmetric case. For simplicity we restrict the 
study of the mass dependence of T c in the first order transition region to the SU(3)- 
symmetric case. In Fig. 7 we see an increase of T c with the current quark mass 
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m [MeV]. For realistic mass values of m w = 129.3 [MeV], m K = 490.7 [MeV], 



Figure 7: The critical temperature T c as a function of the average light quark mass 
m [MeV] in the S"?7(3)-symmetric case. For further explanations see the text. 

m v = 544.7 [MeV] (not depicted in Fig. 7), the rapid crossover sets in at T ~ 181.5 
[MeV] and becomes slow at ~ 192.6 [MeV]. We associate a critical temperature "T c " 
~ 187 ± 0.5 [MeV] to this crossover phenomenon for comparison with extrapolated 
values of T c in chiral perturbation theory (in our case T c is localized as the point of 
inflection in the crossover curve). Thus T c has increased by 5.3 % compared to the 
chiral limit. This result is in agreement with the estimate of Leutwyler P3fl , who 
predicts AT C /T C ~ 5% if realistic pion masses are substituted for the chiral limit 
with 77V = 0. When the critical temperature is extrapolated in the framework of 
chiral perturbation theory, the inclusion of finite quark masses delays the melting 
of (qq) by AT ~ 20 [MeV] ||], while the inclusion of heavier mesons in a dilute 
gas approximation has the opposite effect. At finite quark masses it accelerates the 
melting from T c ~ 240 [MeV] to 190 [MeV]. The delay in the melting due to finite 
quark masses is intimately related to the size of the latent heat. The relation is 
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revealed in the derivation of a Clausius-Clapeyron-equation for QCD, cf. [29]. Thus 
the T c (m)-dependence is conclusive for the strength of a first order chiral transition. 



The case of realistic mass splitting, e^" 11 /e^ lt = 0.77. For the realistic mass 
splitting induced by m s /m U ( i = 18.2 we find for the critical field strengths in large- 
N f 

e° rit < (7 ± 2) • l(T 3 [GeV] 3 , ef* < (-9.09 ± 2.6) • 10~ 3 [GeV] 3 (58) 

with ra % = 68.7, nix = 276.8, m v = 313.4, m v i = 928.6, m CT7r = 887.9, m ai< = 
890.9, m ari = 924.3, m CT , = 698.2 (all masses in units of [MeV]) or an average 
pseudoscalar octet mass of (m ps ) = 203 [MeV]. The critical values for the current 
quark masses are = 2.96 ± 0.85 [MeV] and mf* = 54 ± 15.4 [MeV]. The 

reason why we give an upper bound on e^f rather than precise values for the first 
order transition boundary is the same as in the S*?7(3)-symmetric case. The bound 
on £Q r g* could still be improved by measuring data between e = 4 • 10~ 3 [GeV] 3 , 
Es = —5.2 ■ 10~ 3 [GeV] 3 , where the transition is still of first order, and the above 
values for e^g. It should be remarked that only a rather fine resolution of the 
(qq)/ (qq) (T) curve in steps of AT < 0.1 [MeV] has revealed the true 1st order 
nature of the transition for Eq = 5 ■ 10~ 3 and £$ = —6.5 ■ 10~ 3 [GeV] 3 due to the 
small gap in the condensate. The resolution in units of AT = 1 [MeV] had suggested 
a smooth crossover behaviour already for these smaller values of £o,8- 

Thus we see that also in a large ^/-expansion tiny values for the quark masses 
are sufficient to weaken the chiral transition and turn it finally into a crossover phe- 
nomenon. Similarly tiny current quark masses are sufficient to eliminate the chiral 
transition in an 5£7(3)-Nambu-Jona-Lasinio model as a function of temperature and 
nuclear density |30[] . 

The large- iVy-results for E^f are clearly above the mean-field values of Eq. (f47|). The 
result is plausible as our saddle-point approximation goes beyond the mean-field cal- 
culation. The leading term in a 1/ ^/-expansion corresponds to the summation of a 
class of diagrams, called "super-daisies" B1J . Super-daisies have been summed up by 



Dolan and Jackiw |2B[] to circumvent the IR-divergence problem in an iV-component 
(/> 4 -theory. Our application of the 1/Np- method to the linear sigma-model has been 
similar in spirit. Certainly we cannot claim that the fluctuations we have included 
so far are representative for all fluctuations. In fact, the classical cubic term of our 
potential may still dominate the driving mechanism for the first order transition 
below the critical field strengths. Only the lattice calculation includes all fluctua- 
tions by simulating the full partition function at once (at least in principle). This 
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may explain, why a discrepancy of m 0.25 remains between the large-iV/-ratio and 
the lattice result for m^/m U} d- We should keep in mind, however, that the lattice 
result itself is not yet reliable, as we have argued in section 1. 

The next question which arises in a comparison with the mean-field calculation 
concerns the critical renormalized masses which should vanish for the critical field 
strengths at a second order transition. In the mean-field calculation m? vanishes 

n' 

by construction for the set of critical parameters, and it is easily checked that m a , 
is the only mass that vanishes at criticality. Our value for m a , = 749.5 [MeV] is 
the tree level input mass at zero temperature, which need not vanish. The effective 
masses Xq are temperature dependent. They slowly de/increase with T as we will 
see in section 6. We have mentioned in the discussion of the chiral limit (cf. section 
4) that Xq, Q = 1...4 violate Goldstone's theorem, if they coincide with a dynamical 
mass in our approximation. Similarly it is here not obvious that one of these masses 
should induce an infinite correlation length as T approaches T c . A careful renormal- 
ization prescription should be imposed to identify the renormalized mass(es) that 
go to zero for critical external fields. It would further allow an identification of the 
universality class of the SU(3) x SU( 3) -linear sigma-model for critical parameters 
//q, fi, f2, g, Bo, Eq. We will investigate these questions in a forthcoming work. 

Critical meson masses for m s = 0. Here the large- iV/-results are even somewhat 
smaller than the mean- field results (cf. Eqs. (|55|)). We find 

£ cnt = ( 4 ± !) . i(T 4 [GeV] 3 £% u = (3.06 ± 0.76) ■ 10- 4 [GeV] 3 . (59) 

The associated average pseudoscalar octet mass is (m ps ) = 57.2 ± 28.8 [MeV], and 
m u,d = 1-7 ± 0.4 [MeV], while m s = by construction. 

The Columbia plot (cf. Fig.l) has suggested a concave shape of the first order 
transition boundary. The three critical masses in mean-field are compatible with 
such a shape in an (m^m^)- or an (m, m s ) -diagram, although one should keep in 
mind the large error bars due to the ambiguity in the tree-level parametrization of 
the sigma-model and the sensitive dependence of the boundary on the o^^-mass 
input. The three critical masses in large- Nf do not confirm the conjectural concave 
shape. Only for a realistic ratio of Eq/es the large- iVj-result lies clearly above the 
mean-field value for the critical masses. The error here is at least as large as in the 
mean-field case. The error could have been further reduced, but the remaining size 
would reflect the difficulty in disentangling a very weak first order transition from 
a rapid crossover phenomenon as we have mentioned above. 

In the next section we will see the change in the crossover behaviour, as Eq,Es 
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are further increased to induce realistic mass values. 



6 The realistic mass point 

If we use for /i 2 ,, fx, f 2 , g the values of the chiral limit and choose Eq = 0.0265 [GeV] 3 , 
Eg = —0.0345 [GeV] 3 , Eqs. ( |Tti| - |TT| ) lead to (pseudo) scalar meson masses which are 
listed in Table 1. A comparison to the experimental values shows reasonable agree- 
ment for the pseudoscalar mesons. Therefore we call this point the "realistic" mass 
point. The experimental values which are associated to the scalar meson masses de- 
pend on the identification, which is indicated in a separate row of Table 1. The mass 
splitting between Kq and ao comes out too small in our case. We could have further 
optimized our choice of Eq and e$ to improve the agreement with the experimental 
mass values, but such an optimization should be inconsequential for our results. 



Table 1: Tree level parametrization of the SU(3) x SU(3) linear sigma- model for 
the realistic mass point 

Input 



Aig [GeV] 2 


h 


h 


9 [MeV] 


U [MeV] 


e [GeV] 3 


e 8 [GeV] 3 


5.96- 10~ 2 


4.17 


4.48 


-1,812.0 


94 


0.0265 


-0.0345 



Output (all masses are understood in units of [MeV]) 





m n 


m K 








m* K 






real. 

mass 

point 


129.3 


490.7 


544.7 


1045.5 


1011.6 


1031.2 


1198.0 


749.5 


exp. 

mass 

values 


138.0 


495.7 


547.5 


957.8 


980 if 
'iii 


1322.0 if 
a K = 
K* 


1476.0 if 
a v = 
/o(1476) 


975 if 

Cty = 

/o(975) 



Crossover in the condensates. The crossover behaviour for the normalized light 
and strange quark condensates are displayed in Fig. 8. The rapid crossover leads 
to a decrease of ~ 50 % of the condensate at zero temperature (qq)o in (qq)r over 
a temperature interval AT ~ 10 [MeV], while (ss)t stays remarkably constant up 
to T ~ 197 [MeV] where it starts to decrease rather slowly. The physical reason is 
obvious. It is harder to excite mesons with strange quarks than with light quarks. 
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Figure 8: Light ((qq)) and strange ((ss)) quark condensates normalized to their cor- 
responding values at zero temperature as a function of temperature. The crossover 
behaviour is most rapid between 181.5 < T < 192.6 [MeV]. 



The same qualitative behaviour of (ss)t has been noticed by Hatsuda and Kunihiro 
32[ in a U^ f (3)-version of the Nambu-Jona-Lasinio model. Also the location of the 



crossover region in the NJL-model is around T ~ 200 [MeV]. 

We have indicated the error bars only in the crossover region where they are largest. 
As outlined in section 3, their main sources are the uncertainty in precisely locating 
the saddle point value sad* and the minima (o" ), (cr 8 ) of U e ff- When the errors 
in the current quark masses are added, which enter Eq. (14), we obtain A(qq) ~ 
3.1 ■ 10- 3 [GeV] 3 ,A(ss) ~ 6.7- 10~ 3 [GeV] 3 at T = 180 [MeV] in the transition 
region. In Fig. 8 we have indicated only the numerical errors, the contribution from 
the current quark masses has been left out. Compared to critical meson masses 
(eg™* = 7- 10~ 3 [GeV] 3 ,£g™* = —9.1 ■ 10~ 3 [GeV] 3 ), the crossover happens over a 
larger temperature interval. We find A[(gg) T /(gg) ] rea YA[(gg) T /(gg) ] c "* ~ 52 % if 
^[{qq)r I {qq)o\ denotes the normalized condensate change per 1 [MeV] temperature 
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interval in the rapid part of the crossover region. Nevertheless the crossover in the 
quark condensate (qq) seems to be sharp even for realistic masses. Such a rapid 
change may lead to visible changes in hadron masses which depend on temperature 
and condensates. It could be manifest in hadronic or dilepton spectra in heavy-ion 
experiments. 



Figure 9: Contour plot of U e ff(cr ,as) in a projection on the (a , cr 8 )-plane for four 
temperatures. In Fig. 9a) (T = 10 [MeV]) we see two separated minima (indicated 
with crosses), one unphysical with a < 0,a 8 < and the absolute physcial mini- 
mum with o"o > 0, o"8 < 0. In the crossover region at T = 190 [MeV] (Fig. 9b)) 
both minima lie in the physical sector with <7o > 0, <j$ < 0. The minima coincide for 
temperatures > 200 [MeV] (Fig. 9c) with T = 200 [MeV]). For increasing temper- 
ature the one remaining minimum moves inwards, but stays away from zero due to 
£ ,8 ^ 0, as is seen in Fig. 9d) for T = 350 [MeV]. 

The crossover behaviour in the condensates is also revealed in two-dimensional 
contour plots of the effective potential U e ff(a , a 8 ; sad*). In Fig. 9 we show equipo- 
tential lines for four temperatures in a projection of U e ff on the (<jq, o"8)-plane. In 
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Fig. 9a) we see two separated local minima for T = 10 [MeV], one in the sector 
do < 0, as < and another one in the physical sector o~o > 0, erg < 0. (In the chiral 
limit we would find three minima due to Eq = = In the crossover region 
at T = 190 [MeV] (Fig. 9b)) we still have two separated minima, where the left 
one is a local minimum and the right one is the absolute minimum of U e ff- Now 
both minima lie in the physical sector. At T = 200 [MeV] (Fig. 9c)) the minima 
have merged in a single one. As the temperature is further increased, the absolute 
minimum moves inwards, but stays away from zero due to the explicit symmetry 
breaking fields Sq and e$ (cf. Fig. 9d) for T = 350 [MeV]). The elliptical shape of 
the contour plot reflects the violation of spherical symmetry in (cr , cr 8 )-space due 
to the explicit symmetry breaking, the cubic term in the classical potential and the 
thermal part. As the thermal part gets more dominant at higher temperatures, its 
high-temperature limit changes the elliptical shape to an approximately spherical 
symmetric shape. The two minima below 200 [MeV] may be easily misinterpreted 



as the coexistence of the symmetric and the chiral phase |18[], if one does not realize 
that the minimum (<Jo, as) on the r.h.s. of Fig. 9b) stays the absolute minimum for 
all temperatures. 

Thermodynamics. Further characteristics of the crossover phenomenon are the 
variations in the energy and entropy densities over the temperature interval in the 
crossover region. At this place one should recall the very definition of a first order 
phase transition. At the transition point at least one of the first derivatives of 
a suitable thermodynamic potential should behave discontinuously in the infinite 
volume limit. Thus a crossover in the condensates in general does not exclude a 
finite gap in the energy or entropy densities. In Fig. 10 we have plotted s/T 3 ,e/T 4 , 
and p/T 4 as calculated according to Eqs. ([Z7|-|2"9l) . The data points are strongly 
fluctuating within a range which is indicated by the error bars for the entropy 
curve. The errors for the energy density are of similar size. Errors enter via the 
effective masses Xq, which depend on (cr ), (cr 8 ) and sad*. We have used for A(o~ ) ~ 
(1.0 - 3.0) • 10~ 3 [GeV]= A(cr 8 ) and A(sad*} ~ 0.11 [GeV] 2 in the crossover region 
between T = 182 — 193 [MeV]. The pressure behaves continuously as a function of 
T if it is calculated as p = (—U e ff)- There is only a change of slope in the critical 
crossover region. A direct calculation of p with an integral expression pretends a 
discontinuity. The p/T 4 -curve in Fig. 10 is obtained from (p = — U e ff)- As change 
in the entropy density we find from the actually measured data 

T ■ As < 0.16 [GeV/fm 3 ] (60) 

in a temperature interval 181.5 < T < 192.6 [MeV], where T ■ As is calculated as 
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Figure 10: Entropy density s over T 3 , energy density e over T 4 and pressure p over 
T 4 in the large iVy-expansion for the realistic mass point. Errors are only indicated 
for s/T 3 . 



T 2 • s(T 2 ) — Ti • s(Ti). As rapid change in the energy density we find 

Ae < 0.13 [GeV/fm 3 ] (61) 

or Ae/T 4 = e(T 2 = 192.6 [MeV])/T 4 - e{T x = 181.5 [MeV])/T 4 = 0.29 over the 
same temperature range. For comparison we mention that the gap in the gluonic 



energy density Ae g in a pure SU (3) gauge theory leads to 0, |34}1 

Ae g _ f 2.44 ± 0.24, for AT T = 4 
"T^ - I 1.80 ±0.18, for A^ T = 6, 



(62) 



where A^ T refers to the lattice extension in time direction. These values are by an 
order of magnitude larger than our value for the mesonic contribution Ae/T 4 ~ 0.25, 
defined as indicated above. The decrease of Ae g /T^ under an increase of the time 
extension A^ T indicates strong finite size effects. Going to larger lattices this tendency 
may continue and further reduce the latent heat, but it also gives us a warning. The 
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contribution of Ae g to the total energy gap may be superimposed on the slow change 
of e that we have found in the crossover region and make the crossover in the total 
energy density more rapid. If the size of the errors in sad*, (cr ), (a 8 ) is assumed as 
above, s could vary in the crossover region between T\ and T 2 like 



s(T 2 ) • T 2 - s{T x ) ■ T x < 0.18 [GeV/fm 3 



(63) 



where the resulting error A(s/T 3 ) has been estimated as ±0.46, cf. Fig. 10. In 
an infinitesimally small temperature interval such a gap in s would lead to a finite 
latent heat of < 0.2 [GeV/fm 3 ]. Thus Eq. ( |63| ) gives a loose upper bound on the 
latent heat which could be compatible with our data due to the large errors in the 
crossover region. The bound comes out even smaller if the error of 0.46 is interpreted 
as As/T 3 with "T c " = 187 [MeV]. It leads to AL = T c - As < 0.074 [GeV/fm 3 ]. Both 
bounds are even smaller than Leutwyler's value of 0.4 [GeV/fm 3 ] for T • As [29|, 
which has been obtained from a Clausius-Clapeyron equation in the framework of 
chiral perturbation theory. The small size of the latent heat is finally a consequence 
of the sensitivity of T c to the inclusion of finite quark masses. 
In view of heavy-ion experiments there need not be a latent heat in the strict sense, 
which occurs over a time period during the phase conversion, where the temper- 
ature stays exactly constant. Multiplicity fluctuations in rapidity distributions of 
charged particles are also induced if As is sufficiently large over a small, but finite 
temperature interval (say of the order of ~ 10 [MeV]). In van Hove's formulation 
IH] the physical conditions of a first order transition are identical to those of a rapid 



crossover phenomenon. 

Clausius-Clapeyron equations relate the discontinuities in the condensate and 
the entropy/energy densities. Although they strictly apply to first order transitions 
in the form of Eqs. 



)4|-^) below p9[ , we have tested the relations for our rapid 
crossover phenomenon in two forms. The first one is 

AT r 



disc 



Am 



disc s 



(64) 



where disc...\ r , c , refers to the rapid change ("discontinuity") in the crossover region 
and AT C is the change in "T c " under a variation Am of the current light quark 
masses. The second version is 



(65) 



AT C 




disc{qq)T 






T c 


r.c. 




r.c. 


Ae 



In Eq. (H) we use AT J Am = (0.187 - 0.178)/(0.011), As = 0.0061 [GeV] 3 , 
disc(qq)T = 0.5 ■ (0.22) 3 , and obtain 0.005 for the l.h.s. and for the r.h.s. of Eq. (164 
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In Eq. ( |65"D we have AT C as above, T c = 178 [MeV], disc{qq)p/ (qq)o = (9<z)ti/ (<?<?)o — 
(QQ)t 2 / (qq)o — 0.5, f n = 94 [MeV], m, = 129 [MeV] (the value of our realistic mass 
point), and Ae = 0.001 [GeV] 4 . This way we obtain for the r.h.s. 7 % and for 
the l.h.s. 5 % in Eq. (|65|). The agreement of the order of magnitude on both sides 
indicates that the crossover phenomenon is still rapid enough to satisfy analogous 
relations to a first order transition in the strict sense. 

The same relations have been checked for the first order transition in the chiral 
limit. Eq.s (64), (65) with the appropriate chiral input data predict for Ae/T c 4 ~ 
1.0 ± 0.2 ~ As/T c 3 . The actually measured gap As/T c 3 in chiral thermodynamics 
comes out as 0.5 ±0.2. An error of 0.2 is easily induced by an error of 10 [MeV] in the 
condensates (<7 ), (cr 8 ) at the phase transition. Thus the relations are approximately 
satisfied also in the chiral limit within the relatively large errors. 

We have further analyzed the contribution of the strange (pseudo) scalar mesons 
K and ok to the total energy density e. Their contribution can be completely 
neglected below 40 [MeV]. It increases with temperature to ~ 31 % in the crossover 
region around T = 187 [MeV]. After the crossover the tendency continues, but is 
no longer conclusive for us due to the lack of quark degrees of freedom in the chiral 
symmetric phase. 



The strangeness content of the plasma has been estimated in a lattice simulation |35 
with light quark masses of m u /T = rria/T = 0.05 and one heavier quark mass of 
m s /T = 1.0. In the transition region one finds for the ratio of the fermionic energy 
densities Sp(m s /T = l)/ep(m u /T = 0.05) ~ 0.5. The good agreement with our 
ratio of mesonic contributions from strange and non-strange (pseudo) scalar mesons 
(^(strange mesons) /e(non- strange mesons))~ 0.45 may be accidental, because the 
lattice estimate is based on perturbative relations for the energy density. A fully non- 
perturbative lattice calculation along the lines of Engels et al. |33| in a pure SU(3) 



gauge theory is still outstanding if fermions are included. Nevertheless, it should be 
mentioned that an extrapolation for fermionic energy densities epimjT = 1) and 
Ep{mjT = 0) from the plasma phase to the transition region under the assumption 
of ideal fermion gases leads to Sp(m/T = 1) / 'epimjT = 0) = 0.88 []3"6|1 , which is 
twice the amount of our ratio for mesons in the transition region. This may be taken 
as an indication that the strangeness contribution to e around T c cannot be derived 
from an underlying ideal gas of quarks of m s and m Ui rf-quarks. 

Temperature dependence of effective masses. Fig. 11 displays the temperature 
dependence of the effective masses Xq, Q — 1, 8 up to temperatures of 300 [MeV]. 
We have found a degeneracy of all (pseudo) scalar meson masses at higher tempera- 
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Figure 11: Temperature dependence of the effective mass squares Xq, Q = 1, 8 for 
realistic meson masses. From bottom to top: X n , X K , X v , X a , X Cv , X CK , X v i,X ari . 

tures (T > 400 [MeV]), where the sigma- model fails to describe the plasma phase. 
Immediately above the rapid crossover region chiral symmetry is partially restored. 
The restoration in strange and non-strange sectors is achieved quite differently. We 
find an approximate degeneracy between X n and X a ,, Xx,X v and X aK ,X r] r and 
X a7T , while X a becomes degenerate with other masses only at T ~ 400 [MeV] . The 
ovmass never meets the 7r-mass, because the [Tritanomaly term (proportional to 
g) enters the mass formulas with opposite sign. From lattice results it is very likely 
that other modes than X n , X a , dissolve in their constituents and cease to exist in 
the chiral symmetric phase. Coming from the broken phase we notice that the effec- 
tive mass X a , drops to a remarkably small value in the crossover region {X c , ~ (X n 
at T = 0) ~ 130 [MeV]), while X^ is almost zero. Furthermore, we have displayed 
the errors in X n and X a , in the transition region. The errors have increased with T 
to these values due to the contribution of Asad(T). For X n the error is particularly 
large in the crossover region, since the error AX n is proportional to l/X w and X n is 
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almost zero in this region. Before we jump to conclusions about measurable effects 
due to tiny a v ' and 7r-masses in the crossover region, we should clarify the precise 
physical meaning of the effective masses Xq. We recall that the main contribution 
to the temperature dependence of Xq comes from the saddle point contribution 
sad*(T), which is just the leading term in the 1/iVj-expansion. 

7 Summary of results and conclusions 

In agreement with the general expectation results of the high-T-expansion are qual- 
itatively correct, but fail quantitatively. The critical temperature in the chiral limit 
deviates by ~ 85 [MeV] from T c in the numerical evaluation, which is applicable also 
in the transition/crossover region. The crossover region for realistic meson masses 
is shifted by roughly 80 [MeV] between the high-T- and the numerical results. In 
our model the high-T-expansion practically never reproduces the numerical results 
in a quantitative way. While the meson condensate (<7 ) in the realistic mass case 
has dropped to values < 2 [MeV] around T ~ 400 [MeV] in the high-T-expansion, it 
is still larger than 30 [MeV] at this temperature in the numerical calculation. The 
temperature lies already far outside the applicability range of the model. 

The crossover in the light quark condensate looks still rapid for realistic meson 
masses (A((qq) T ) ~ 50% of (qq)o in a temperature interval of 10 [MeV] in the 
crossover region). Less rapid look the variations in the energy and entropy densities 
between 181.5 [MeV] < T < 192.6 [MeV]. We find Ae ~ 0.13 ± 0.02 [GeV/fm 3 ] 
and T ■ A s ~ 0.16 ± 0.02 [GeV/fm 3 ]. As a loose bound on the latent heat we 
obtain 0.2 [GeV/fm 3 ], as a more stringent one 0.1 [GeV/fm 3 ], if the change in 
e and s would occur over an arbitrarily small temperature interval. The main 
contributions to the numerical errors which prevent us from a unique interpretation 
of the transition/crossover region come from the uncertainties in the saddle point 
value and the minima ((Jo), (ct 8 ) of U e ff- 

For temperatures above 116 [MeV] the saddle point variable gives a positive 
contribution to the effective masses entering the argument of the logarithm in U e ff. 
It increases with temperature. This was one of the reasons why we have chosen the 
large- iVj-expansion. The original hope to completely avoid imaginary parts of U e ff 
in this scheme could not be confirmed by the results. The effective potential is still 
plagued with imaginary parts for certain regions, where \<Jq$\ < | (co,s) I an d over the 
entire temperature range we have studied (up to 250 [MeV] and above). 
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The large error bars on e and s in the transition region may leave some doubts on 
the smooth nature of the conversion between the chiral symmetric and the chirally 
broken phase, in particular, because a smooth crossover in the condensates does not 
automatically imply a smooth change in e and s. Even if the transition would be of 
first order, and even if we use a loose bound on the latent heat, it is as small as 0.2 
[GeV/fm 31 



The ratio of critical to realistic current light quark masses has been estimated 
as m^rf/m^d ~ 0.03 ± 0.02 in mean-field and as / Tn u ,d ~ 0.26 ± 0.08 in large- 
Nf. The large- Nf result implies that realistic quark/meson masses lie less deeply 
in the crossover region than the mean-field result suggests. Due to the fluctuations 
which are effectively included in the large- Nf- approximation the ratio in large- Nf 
is about half of the lattice result. The average pseudoscalar octet masses are 126.4 
[MeV] and 203 [MeV] for the realistic ratio of e^ 11 je c { %t in mean-field and large- Nf, 
respectively. Note that an error of 30 % in the critical quark masses blows up to 
an error of 55 % in the critical meson masses, if the meson mass squares are linear 
in the current quark masses. Thus the critical values for the meson masses are not 
very conclusive. 

In view of heavy-ion experiments one may conclude that there is still little hope 
for measurable experimental signatures of the crossover region, although the ratio 
rn ud/ m u,d has become larger than the mean-field estimate. Quantities which are 
indirectly observable such as condensates, energy and entropy densities vary too 
smoothly during the phase conversion. Experimentalists should be warned, however, 
to accept these conclusions without care. 

The tendency that m^a / 'ra u ,d is larger in large- Nf than in mean-field is plausible 
as an effect of the included fluctuations, but as we have seen this effect could not 
be confirmed for the S{7(3)-symmetric or the m s = case. The critical mass ratios 
there were even smaller in large- Nf than in mean-field or at least of the same order 
of magnitude within the errors. This gives a hint that the fluctuations we include 
in our approximation are likely not the only important ones. In particular, it is 
not clear that they account for fluctuations which induce a renormalization of the 
quartic and cubic couplings in the Lagrangian. 



Recently it has been raised by Gavin, Goksch and Pisarski |I7J that the first 



order of the chiral transition may be mainly fluctuation induced. A fluctuation 



induced transition has been first discovered by Coleman and Weinberg |R7[. It 



refers to the situation that a system with more than one relevant coupling can have 
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a first order transition induced by quantum fluctuations. In an understanding of the 
renormalization group the fluctuation induced transition occurs due to so-called run- 



away RG trajectories (see e.g. [p8|). Such an origin for a first order transition has 
been demonstrated in 4 — e dimensions for the SU(N) x SU (N) linear sigma model 
by Paterson f3l| , for more general symmetry groups (O(N) x O(M), U(N) x U(M)) 
by Pisarski and Stein [PC[j . In 4 dimensions Shen |38[] has shown in a numerical 
simulation of a U(N) x U(N) symmetric scalar model that a fluctuation induced 
first order transition occurs in particular for large coupling f'2. For large and slowly 
varying f2, fx runs fast from large to small values as a function of energy scale. The 
same tendency in an SU(3) x SU (3) linear sigma model would be compatible with a 
larger critical mass ratio m c ™J{ / m u ^- The larger mass ratio of the lattice simulation 
points into the same direction, since the lattice includes all fluctuations at once. 
In our present approach it is difficult to disentangle the driving mechanism for the 
first order transition below the critical mass values. A further improvement by 
including subleading corrections in 1/Nf or an e-expansion in d = 4 — e dimensions 
should clarify, whether the ratio m c ™d/m u ^ does change to larger values. If the 
realistic finite quark masses lie in fact closer to the first order transition boundary 
than our results suggest, remnants of the first or second order transition may be 
more easily visible in experiments. 

A second warning should be mentioned not to take the smooth crossover for 
guaranteed by the present approach. So far we have completely neglected the quark- 
and gluonic substructure. In particular the rearrangement of gluonic degrees of 
freedom has not been taken into account. The gluonic contribution to the change 
in entropy and energy densities may well accelerate the crossover process. Note 
that we have chosen the couplings in the sigma-model as temperature and energy 
scale independent over all temperatures up to T c . In principle, the temperature and 
scale dependence of the couplings should be derived from QCD rather than being 
assumed. Less ambitious it may be derived from an effective model underlying the 
sigma-model and containing quark and gluonic degrees of freedom. For the critical 
mass ratio this offers at least a chance for being closer to 1. Work in this direction 
is in progress. 
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